System and method for estimating the multi-path delays in a signal using a spatially blind antenna array

ABSTRACT

A method is described that enables a system to estimate the individual times of arrival of multi-path signals components in a received signal while blind to the particular angular response characteristics of an antenna array. Additionally, a system is described that estimate the impulse response of the modulation channel. The impulse response is used to generate impulse response estimates for all multi-path signal components without reference to the actual angle of arrival of the signal. The impulse response for the multi-path constituents are directly associated with the time delays on each path, providing delay estimates for each path.

CROSS REFERENCES

The present application is with and claims priority benefit of U.S. provisional application entitled “Geolocation of Mobile Appliances” Ser. No. 60/418,342 and filed on Oct. 16, 2002, the entirety of which is incorporated herein by reference.

The present application is related to and concurrently filed with applications titled “A NETWORK OVERLAY GEO-LOCATION SYSTEM WITH SMART ANTENNAS AND METHOD OF OPERATION” Ser. No. 10/351,040, “WIRELESS COMMUNICATION NETWORK MEASUREMENT DATA COLLECTION USING INFRASTRUCTURE OVERLAY-BASED HANDSET LOCATION SYSTEMS” Ser. No. 10/531,042, “NETWORK OVERLAY LOCATION SYSTEM AND METHOD FOR AIR INTERFACE WITH FREQUENCY HOPPING” Ser. No. 10/531,041, “A SYSTEM AND METHOD FOR ENHANCING THE ACCURACY OF A LOCATION ESTIMATE”, Ser. No. 10/531,044, and “SYSTEM AND METHOD FOR OPERATING A NETWORK OVERLAY GEO-LOCATION SYSTEM WITH REPEATERS” Ser. No. 10/531,038, filed Oct. 16, 2003, the entirety of each of these applications is incorporated herein by reference.

BACKGROUND

Digital signals are often filtered using a pulse shaping filter prior to transmission. This is typically done to contain the signal bandwidth and minimize intersymbol interference between signal components corresponding to different digital symbols. This is shown in FIG. 1, where symbols from a digital constellation corresponding to the information being modulated are passed though a pulse shaping filter. In FIG. 1, the digital data to be transmitted is mapped into a complex signal constellation in block 101. For example, the complex signal constellation used may be an M-ary QAM constellation; however other constellations are also used. The mapped constellation undergoes pulse shaping in a filter as shown in block 102. Several methods known in the art can be employed for pulse shaping. The filtered constellation signal is converted to a radio frequency, represented as block 103, for transmission over the ether.

An artifact of this signal generation provides an unintended feature that has been the focus of recent investigation. If the digital signal is passed through a multi-path channel, the channel output is an aggregate of delayed, possibly faded and phase shifted replicas of the original digital signal. In practice, this occurs if a multiplicity of reflections of the transmitted signal are contained in the received signal. The delays can either be absolute, if the time of arrival of the direct path signal is known, or could be relative delays between the multi-path components. If these multi-path signals are received at the antenna array, the received signal can be mathematically formulated as a space-time signal.

When the characteristics of the pulse shaping filter and the antenna array are known, a theory of signal processing can be applied to estimate the multi-path delays of the signal components and their particular directions of arrival. This signal processing analysis is referred to as space-time processing. Space-time processing is a group of techniques that may be applied to resolve the received space-time signals into a sum of faded space-time signals. Each of these space time signals corresponds to the particular angle of arrival and time delay of one of the original multi-path signal components.

It is advantageous to develop a mathematical description of the prior art technique to convey the manner in which multi-path delays and angle of arrival (“AOA”) are currently calculated to fully appreciate the distinctness of the to be disclosed subject matter. The prior art method is illustrated in FIG. 2.

A column vector r_(k) denotes the received signal at antenna k of an antenna array with m antennas, where k=1,2, . . . m. An impulse response h_(k) of the multi-path channel is derived from r_(k), represented in block 201. The derivation of the column vector h_(k) can be achieved by various methods and implemented with signal processors through software and/or hardware.

If the source data associated with this received block is known, a simple means of extracting h_(k) is via the delay matrix corresponding to this source data. The delay matrix Z is formed by stacking symbol shifted copies of the source data in rows to a depth that defines the extent of the desired impulse response and truncating its longer dimension to match the length of r_(k). An estimate of h_(k) is given by: h_(k)=(ZZ^(H))⁻¹Zr_(k).

Alternate means for estimating the impulse response may provide better or worse estimates, depending on the particular modulation format of the data, the block length, the fading characteristics of the multi-path channel and possibly other parameters. Some of these other methods are blind to the actual data transmitted, using properties of either the signal modulation and/or of the channel instead.

Having estimated the impulse response of the multi-path channel from the source to antenna k of the array, a vectorized space-time impulse response over the entire array is formed, in block 202, by stacking the individual impulse response estimate h_(k) into a long column vector {right arrow over (I)}, given by:

$\overset{\rightarrow}{I} = {\begin{bmatrix} h_{1} \\ \ldots \\ h_{k} \\ \ldots \\ h_{m} \end{bmatrix}.}$

Theoretically, {right arrow over (I)} can be expressed as

${I = {{\sum\limits_{i = 1}^{n}I_{i}} + N}},$ where {right arrow over (I)} indexes the individual space-time impulse responses, i=1,2, . . . , n, of the distinct multi-path components and N is a noise vector.

Any particular I_(i) is of the form I_(i)=β_(i)ηa(θ_(i)){circle around (×)}g(τ_(i)) in which β_(i) denotes the fade multiplier for the signal block, and η denotes the signal amplitude at the transmitter. a(θ_(i)) denotes the antenna response corresponding to a signal arriving from angle θ_(i), {circle around (×)} denotes the Khatri-Rao product, and g(τ_(i)) denotes the pulse shaping waveform delayed by τ_(i) and sampled. This formation of an outer product and aggregate in the covariance matrix, is represented in block 203. The formation of the covariance matrix can be implemented with signal processors or other computer processors through software and/or hardware devices.

The vectors a and g can be expressed as:

${{a\left( \theta_{i} \right)} = \begin{bmatrix} {a_{1}\left( \theta_{i} \right)} \\ {a_{2}\left( \theta_{i} \right)} \\ \ldots \\ {a_{m}\left( \theta_{i} \right)} \end{bmatrix}},{and}$ ${g\left( \tau_{i} \right)} = {\begin{bmatrix} {g\left( {{{- l}\; T_{s}} - \tau_{i}} \right)} \\ {g\left( {{{- \left( {l - 1} \right)}T_{s}} - \tau_{i}} \right)} \\ \ldots \\ {g\left( {{l\; T_{s}} - \tau_{i}} \right)} \end{bmatrix}.}$

In the equation for g(τ_(i)), l denotes the sampling depth of the pulse shaping function and T_(s) is the sampling time.

Given this formation of the space-time impulse response, when the number of multi-path components is smaller than the dimension of the symbol sampled impulse response vector I, it is possible to estimate the multi-path delays τ_(i) and the multi-path arrival angles θ_(i).

The prior art approach to estimating the delays and arrival angles relies on an explicit knowledge of the aggregate of all vectors:

${a\left( \theta_{i} \right)} = \begin{bmatrix} {a_{1}\left( \theta_{i} \right)} \\ {a_{2}\left( \theta_{i} \right)} \\ \ldots \\ {a_{m}\left( \theta_{i} \right)} \end{bmatrix}$ for all angles θ_(i). This aggregate is termed the array manifold, A. It is assumed that the pulse shaping function at the transmitter is known at the receiver. Denoting the aggregate of all vectors as

${g\left( \tau_{i} \right)} = \begin{bmatrix} {g\left( {{{- l}\; T_{s}} - \tau_{i}} \right)} \\ {g\left( {{{- \left( {l - 1} \right)}T_{s}} - \tau_{i}} \right)} \\ \ldots \\ {g\left( {{l\; T_{s}} - \tau_{i}} \right)} \end{bmatrix}$ for all values of τ_(i) as the delay manifold, G , then the quantity K=A{circle around (×)}G represents the space-time manifold.

The observation that I_(i) is contained in K leads to a primary objective of space-time processing: searching the manifold K for weighted linear combinations of vectors I_(i) such that a best fit to the observed space-time impulse response {right arrow over (I)} is generated as shown in block 204. A variety of techniques may be applied for this purpose, such as Multiple Signal Classification (MUSIC), The Method of Alternating Projections (APM), etc, which can be implemented through software and/or hardware. Other mathematical descriptions for jointly estimating the angle of arrival (“AOA”) and time delays can be found in Ziskind, I., Wax, M., “Maximum likelihood localization of multiple sources by alternating projection”, IEEE Trans. Acoust., Speech, Signal Process. vol. 36, no. 2 (October 1988), 1553-1560; Van Der Veen, M, Papadias, C. B., Pautraj, A. J., “Joint angle and delay estimation” IEEE Communications Letters vol. 1-1 (January 1997), 12-14; Schmidt, R. O. “Multiple emitter location and signal parameter estimation” Proc. RADC Spectrum Estimation Workshop, (March 1999), 243-258; Young-Fang Chen, Michael D. Zoltowski “Joint Angle and Delay estimation of DS-CDMA communication systems with Application to Reduced Dimension Space-time 2D Rake Receivers”, IEEE Transactions on Signal Processing; Paulraj, A. J., Papadias, C. B., “Space-Time Signal Processingfor Wireless Communications”, IEEE Signal Processing Magazine, vol. 11 (November 1997), 49-83; Paulraj, A. J., Papadias, C. B., “Space-Time Signal Processingfor Wireless Communications: A Survey” Information System Laboratory, Stanford University; and Haardt, Brunner and Nossek “Joint Estimation of 2-D Arrival Angles, Propagation Delays, and Doppler Frequencies in Wireless Communications”; all of which are incorporated herein by reference.

An object of the disclosed subject matter is to obviate the deficiencies of the prior art by removing the dependency of the time delay estimates from the spatial and gain characteristics of an antenna array thus allowing multi-path delay estimates to be obtained for any generic antenna array. This object is achieved by recasting the array manifold in a spatially blind manner so as to be independent of the array characteristics.

It is another object of the disclosed subject matter to present an improved method for estimating the multi-path delays in a signal received at any k array element. The method includes estimating an impulse response at each k antenna, generating a space-time impulse response, and forming a covariance matrix and resolving the covariance matrix with a known antenna array manifold. Additionally, a novel improvement to known methods includes the step of resolving the covariance matrix with a fictitious antenna array manifold.

It is still another object of the disclosed subject matter to present a novel method for estimating the multi-path delays in a signal using a spatially blind antenna array. The method includes generating an impulse response h_(k) for each antenna k in the antenna array and determining a vectorized space-time impulse response I over the antenna array. The method further includes creating a covariance matrix C, a fictitious manifold A_(f), where A_(f) is spatially blind and independent of the array characteristics, and then resolving the covariance matrix C with the fictitious manifold A_(f) to estimate the multi-path delays τ_(i) in a manner independent of the array characteristics.

It is yet another object of the disclosed subject matter to present a method of estimating the multi-path delays of a sequence of j blocks of a signal received at an antenna array comprising k antenna elements independently of the spatial array characteristics of the antenna array. The method includes deriving a channel impulse response estimates h_(j,k) for each block j at each antenna k and determining a vectorized aggregate space-time impulse response I for each block j. The method includes the steps of forming an estimated covariance matrix for the sequence of j blocks, forming an array manifold A_(f) void of spatial information; and then resolving the covariance matrix with the fictitious array manifold A_(f) to determine the multi-path delays τ_(i).

It is also an object of the disclosed subject matter to present a novel system for estimating the multi-path delays in a signal using a spatially blind antenna array. The system includes an antenna array, a means for generating an impulse response h_(k), a means determining a vectorized space-time impulse response I and a means for creating a covariance matrix C. The system also includes a means for creating a fictitious manifold A_(f), wherein A_(f) is spatially blind and independent of the array characteristics; and a means for resolving the covariance matrix C with the fictitious manifold A_(f) to estimate the multi-path delays τ_(i) independent of the array characteristics.

These objects and other advantages of the disclosed subject matter will be readily apparent to one skilled in the art to which the disclosure pertains from a perusal or the claims, the appended drawings, and the following detailed description of the preferred embodiments.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a representative flow chart of a portion of a prior art system and/or method for transmitting a digital signal.

FIG. 2 is a representative flow chart of a prior art receiving system and/or method for space-time processing a received signal.

FIG. 3 is a representative flow chart for spatially blind space-time processing a received signal according to a system/method of the present disclosure.

FIG. 4 is a representative flow chart for estimating delays of multipath components of a received signal according to an embodiment of a system/method of the present disclosure

DETAILED DESCRIPTION

The present disclosure presents a system and method for a technique for determining the individual times of arrival of multipath signal components in a received signal. The delays in the multipath components are determined by creating a fictitious array manifold A_(f) in which the directional knowledge of the array is absent. As discussed earlier, the vector a is produced with the knowledge of all the vectors for all angles θ_(i), is represented as:

${a\left( \theta_{i} \right)} = {\begin{bmatrix} {a_{1}\left( \theta_{i} \right)} \\ {a_{2}\left( \theta_{i} \right)} \\ \ldots \\ {a_{m}\left( \theta_{i} \right)} \end{bmatrix}.}$

To remove the dependence on θ_(i), each dimension of the array takes all possible complex numbered values. The independent aggregate vector is then given by:

${{a\left( \theta_{i} \right)} = \begin{bmatrix} a_{1} \\ a_{2} \\ \ldots \\ a_{m} \end{bmatrix}},$ where each a_(k) is contained in the complex number space C, so that the fictitious array manifold A_(f) is the space C^(m). This formulation allows space-time processing to proceed blind to the specific characteristics of the antenna array such as element spacing, element gain, etc.

As discussed above, the signal processing flow for the established prior art technique is shown in FIG. 2. This technique requires knowledge of the antenna array manifold A, i.e., one must have explicit knowledge of the antenna response corresponding to a signal arriving at each possible angle for each element in the antenna array. As discussed earlier, it is not always possible to easily obtain the antenna array manifold A as calibration of the antenna array is a tedious, time-consuming, and expensive task. The prior art technique allows for determination of an estimation of multipath delays and the particular angles of arrival of the multipath signals.

The signal processing flow for the above-discussed subject matter, on the other hand, employs a fictitious array manifold as shown in FIG. 3. An estimate of the impulse response is determined at each antenna element as shown in block 301. From the impulse response, a space-time impulse response is generated in block 302 and an outer product and aggregate in the covariance matrix is formed in block 303 in a manner similar to that as discussed for the prior art. In block 304 the covariance matrix is resolved with the generalized fictitious antenna array manifold A_(f) defined above rather than the known antenna array manifold A as shown in FIG. 2. While the novel techniques described herein does not allow for resolution of the particular angles of arrival of the multipath signals, it does allow for determining an estimation of the time delays of the multipath components given an antenna array with arbitrary and unknown properties. Therefore, the novel techniques herein described allow for the use of antennas separated by arbitrary distances and with arbitrary gain characteristics while still obtaining important information: the time delays of multipath signal components.

To facilitate understanding of the novel techniques described herein, an example follows for the case of a signal received at a two antenna array in which the established prior art space-time techniques derive delay estimates on a multi-path signal, followed by estimates developed by use of a spatially blind antenna array according to an embodiment of the present disclosure.

A sequence of blocks of signal received concurrent in time at the two antennas are denoted by r_(jk) where k=1, 2 and j=1,2, . . . J, the index j counts the blocks and of course the index k references the antennas, two in this case. It is, of course, to be understood by those of skill in the art that the present discussion is exemplary only and that the techniques described herein are in no way limited to antenna arrays composed of only two antenna elements. The disclosed techniques are applicable to antenna arrays of an arbitrary number of antenna elements. If the source data associated with each block is known, possibly from demodulating the received signal at either or both of the antennas, or by demodulation at some other receiver, the delay matrix can be developed as illustrated in the previous section and this can be used to derive channel impulse response estimates for each block at each antenna. Thus h_(j,k)=(Z_(j)Z_(j) ^(H))⁻¹Z_(j)r_(j,k), where h_(jk) is the impulse response estimate for block j and antenna k, and Z_(j) is the delay matrix for block j.

The vectorized aggregate space-time impulse response for block j is given by

${I_{j} = \begin{bmatrix} h_{j,1} \\ h_{j,2} \end{bmatrix}},$ from which an estimated covariance matrix for the sequence of blocks can then be formed as

$C = {\sum\limits_{j = 1}^{J}{I_{j}{I_{j}^{H}.}}}$

Most known methods for extracting an arrival angle do so by operations on C. Other methods that operate on the sequence {I_(j)} differently can also be used.

All established techniques for estimating the path delays and arrival angles in the space-time context require knowledge of the array manifold A which is the aggregate of all possible vectors:

${a(\theta)} = \begin{bmatrix} {a_{1}(\theta)} \\ {a_{2}(\theta)} \end{bmatrix}$ in which a_(k)(θ) is the complex gain of antenna k for a signal arriving from angle θ. Explicit knowledge of antenna gain for each antenna is required. In addition, the pulse shaping function and thus the delay manifold G as detailed in the previous section are also needed.

Most established prior art techniques estimate the delays and arrival angles of the multi-path signals by decomposition of the covariance matrix into a sum of estimated impulse responses. This is done using the space-time manifold: K=A{circle around (×)}G.

The estimated covariance matrix is resolved into estimates of the impulse responses associated with each multi-path component. Each estimated impulse response is drawn from K and therefore points directly to a particular delay and arrival angle.

In accordance with an embodiment of the disclosed subject matter, an example using a spatially blind estimation techniques for the time delay estimates is illustrated below for the same scenarios as described above.

To obtain the delay estimate without knowledge of the antenna array or an explicit knowledge of the antenna gains, a fictitious array manifold is used, A′_(f), given by the aggregate of all vectors,

$a = \begin{bmatrix} a_{1} \\ a_{2} \end{bmatrix}$ where a₁ and a₂ range over the set of complex numbers. It is worthy to note that while the delay estimates can be determined for multi-path signals independently of the antenna characteristics, angle of arrival cannot be determined in the same manner. This is not a problem if all that is of concern is the relative time delays of the multipath signals and the angles of arrival of those multipath signals is not of consequence.

The process estimates the delays without knowledge of the actual array manifold A by using instead the fictitious manifold A′_(f), so that A′_(f) replaces A in the space-time manifold K=A{circle around (×)}G. This then becomes the aggregate of vectors,

${K = \begin{bmatrix} {a_{1}{g(\tau)}} \\ {a_{2}{g(\tau)}} \end{bmatrix}},$ as a₁ and a₂ cover the space of complex numbers and τ covers the expected range of multi-path delays.

The estimated covariance matrix is resolved in order to estimate the impulse responses associated with each multi-path component. Each estimated impulse response is drawn from K_(f) and maps directly to a particular τ and complex vector [a₁, a₂]^(t). The signal processing flow for this embodiment is shown in FIG. 4. An estimate of the impulse response is determined at each antenna element as shown in block 401. From the impulse response, a space-time impulse response in generated in block 402 and an outer product and aggregate in the covariance matrix is formed in block 403. In block 404, the covariance matrix is resolved with the specific fictitious antenna array manifold A′_(f).

The earliest arrival time is of particular interest in many applications such as cellular geo-location. This technique, termed BLind Space Time (“BLST”) permits the estimation of relative delays in multi-path signal components using arbitrary configurations of antennas in an array. This, for example, allows one to apply the method to antennas that are separated by arbitrary distances and with arbitrary gain characteristics. Existing antenna configurations in the cellular wireless network can be used without the need for calibration of the array. Calibration is a tedious task that needs to be updated often as a known array is required for signal processing according to the prior art. Thus the use of Blind-Space Time for estimating multi-path delays allows for the application to any arbitrary antenna array without requiring calibration.

While preferred embodiments of the present inventive system and method have been described, it is to be understood that the embodiments described are illustrative only and that the scope of the embodiments of the present inventive system and method is to be defined solely by the appended claims when accorded a full range of equivalence, many variations and modifications naturally occurring to those of skill in the art from a perusal hereof. 

1. A method for estimating the multi-path delays τ_(i) in a signal using a spatially blind antenna array comprising k arbitrary antenna elements, comprising the steps of: generating an impulse response h_(k) for each antenna element k in the antenna array; determining a vectorized space-time impulse response I over the antenna array; creating a covariance matrix C; creating a fictitious array manifold A_(f), wherein A_(f) is spatially blind and independent of the array characteristics; and resolving the covariance matrix C with the fictitious manifold A_(f) to thereby estimate the multi-path delays τ_(i) independent of the array characteristics.
 2. The method of claim 1 wherein the impulse response estimate h_(k) is determined from the equation: h _(k)=(ZZ^(H))⁻¹Zr_(k) where Z is a delay matrix and r_(k) is the column vector of the received signal at antenna element k of the antenna array, where k=1,2, . . . m.
 3. The method of claim 2 wherein the space-time impulse response vector I is formed by stacking the individual impulse response estimates h_(k) into a column vector.
 4. The method of claim 1 wherein the fictitious manifold A_(f) is the aggregate of all vectors: ${a = \begin{bmatrix} a_{1} \\ a_{2} \\ \vdots \\ a_{m} \end{bmatrix}},$ where a_(k) (k=1,2, . . . m) range over the set of complex numbers, where m is the number of antenna elements in the array.
 5. The method of claim 1 wherein the covariance matrix C is generated according to the following equation: C = ∑I I^(H).
 6. The method of claim 1, wherein the fictitious array manifold A_(f) is used to form the space-time manifold and the space-time manifold operates to resolve the multi-path delays.
 7. The method of claim 1 wherein the step of resolving the covariance matrix C to determine multi-path delays τ_(i) uses the method of MUltiple SIgnal Classification (MUSIC) techniques.
 8. The method of claim 1 wherein the step of resolving the covariance matrix C to determine multi-path delays τ_(i) uses the Method of Alternating Projection (APM).
 9. A method of estimating the multi-path delays τ_(i) of a sequence of j blocks of a signal received at an antenna array of k isotropic antenna elements, independently of the spatial array characteristics of the antenna array, comprising the steps of: deriving channel impulse response estimates h_(j,k) for each block j at each antenna k; determining a vectorized aggregate space-time impulse response I for each block j; forming an estimated covariance matrix for the sequence of j blocks; providing an array manifold A_(f) void of spatial information; and, resolving the covariance matrix with the array manifold A_(f) to determine the multi-path delays τ_(i).
 10. The method of claim 9, wherein the impulse response estimate h_(j,k) for block j is determined from the equation: h_(j,k)=(Z_(j)Z_(j) ^(H))⁻¹Z_(j)r_(j,k) where Z_(j) is a delay matrix for block j and r_(jk) is the column vector of the received signal for block j at antenna k of the antenna array, where k=1,2, . . . m.
 11. The method of claim 10, wherein the space-time impulse response vector I is formed by stacking the individual impulse response estimates h_(jk) into a column vector.
 12. The method of claim 9, wherein the fictitious manifold A_(f) is the aggregate of all vectors: ${a = \begin{bmatrix} a_{1} \\ a_{2} \\ \vdots \\ a_{m} \end{bmatrix}},$ where a_(k) (k=1,2, . . . m) range over the set of complex numbers, where m is the number of an antenna element in the array.
 13. The method of claim 9, wherein the covariance matrix C is generated according to the following equation: $C = {\sum\limits_{j = 1}^{J}{I_{j}\;{I_{j}^{H}.}}}$
 14. The method of claim 9, wherein the fictitious array manifold A_(f) is used to form the space-time manifold and the space-time manifold operates to resolve the multi-path delays.
 15. The method of claim 9, wherein the step of resolving the covariance matrix C to determine multi-path delays τ_(i) uses multiple signal classification techniques.
 16. The method of claim 9, wherein the step of resolving the covariance matrix C to determine multi-path delays τ_(i) uses Alternating Projection.
 17. A system for estimating the multi-path delays τ_(i) in a signal using a spatially blind antenna array comprising: an antenna array for receiving the signal; a means for generating an impulse response h_(k) for each antenna k in the antenna array; a means determining a vectorized space-time impulse response I over the antenna array; a means for creating a covariance matrix C a means for creating a fictitious manifold A_(f), wherein A_(f) is spatially blind and independent of the array characteristics; and, a means for resolving the covariance matrix C with the fictitious manifold A_(f) to estimate the multi-path delays τ_(i) independent of the array characteristics.
 18. The system of claim 17, wherein the fictitious array manifold A_(f) in part forms the space-time manifold and a space-time manifold operates to resolve the multi-path delays. 